# -*- coding: utf-8; mode: snippet -*-
# name: ANOVA/ANCOVA using statsmodels.anova_lm
# key: anova_lm@statsmodels
# contributor: Shuguang Sun
# --
# import statsmodels.api as sm
# from statsmodels.formula.api import ols

${1:model} = '${2:height ~ C(type)}'
${3:anova} = sm.stats.anova_lm(ols($1, data=${4:data}).fit(), type=2)